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We use a simple model of disc-galaxy formation in the ACDM cosmology to simultaneously 
reproduce the slopes, zero-points, scatter and uncorrelated residuals in the observed velocity- 
luminosity (VL) and radius-luminosity (RL) relations. Observed I-band luminosities are con- 
verted to stellar masses using the IMF-dependent relation between stellar mass-to-light ratio, T, 
and color. The model treats halo concentration, spin parameter, and disc mass fraction as inde- 
pendent log-normal random variables. Our main conclusion is that the VL and RL zero-points 
and the uncorrelated residuals can only be reproduced simultaneously if adiabatic contraction is 
avoided. One or the other could be fixed by appealing to unrealistic values for T, halo concentra- 
tion c, and spin parameter X, but not all together. The small VL scatter is naturally determined by 
the predicted scatter in c and in T, quite independent of the large scatter in X. However, the RL 
scatter, driven by the scatter in X, can be as low as observed only if G\ a x ~ 0.25, about half the 
value predicted for CDM haloes. This may indicate that discs form in a special subset of haloes. 
The VL slope is reproduced once star formation occurs only above a threshold surface density, 
and the RL slope implies that the disc mass fraction is increasing with halo mass, consistent with 
feedback effects. A model that incorporates the above ingredients provides a simultaneous fit to 
all the observed features. In particular, the elimination of adiabatic contraction allows 80% disc 
contribution to the observed rotation velocity of bright discs at 2.2 disc scale lengths. The lack 
of halo contraction may indicate that disc formation is not as smooth as typically envisioned, but 
instead may involve clumpy, cold streams. 
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1. Introduction 

The correlations between the global structure parameters of disc galaxies, rotation velocity V, 
exponential radius R, and luminosity L, provide invaluable constraints on the formation process of 
these galaxies. Although various studies have attempted to explain the origin of the RL or the VL 
(also known as Tully-Fisher) relation, the true challenge lies in finding a self-consistent model of 
disc formation that can match both simultaneously. Here we investigate how standard models of 
disc formation fair in this respect. 

We use the observed structural parameters from an /-band sample of ~ 1800 bright disc galax- 
ies compiled by Courteau et al. (2003). In order to facilitate a simple comparison with theory, 
we convert the observed /-band luminosity into a stellar masses M s , using the stellar mass-to-light 
ratio Ti at a given color as obtained by Bell et al. (2003) and Bell & de Jong (2001). The main 
uncertainty in this conversion is the normalization of Tj, reflecting the unknown stellar initial mass 
function (IMF). We adopt the IMF advocated by Bell et al. (2003) and note that the most extreme 
top-heavy IMF realistically considered (the Kroupa-IMF; Kroupa, Tout & Gilmore 1993) has val- 
ues of log(T/) that are only 0.15 dex smaller. We use as joint constraints on our model the observed 
log slopes and zero-points of the conditional V\M S and R\M S scaling laws, the conditional scatter 
about each at a given M, and especially the cross-correlation between the corresponding residuals 
AV and AR at M s . Throughout we assume Ho = 70kms _1 Mpc -1 . 

2. Disc Formation Models 

Our model starts in the spirit of the classical model by Fall & Efstathiou (1980) and Mo, 
Mao, & White (1998, hereafter MMW). Each model galaxy is specified by four parameters: the 
total virial mass M v j r , the halo concentration parameter c, the baryon spin parameter X, and the 
mass fraction that ends up in the disc m&. The dark-matter is assumed to follow a spherical NFW 
density profile, with c that varies systematically with M vu . The disc, with mass M& = m^M Ya , is 
assumed to have an exponential density profile, whose scale length R& follows from X, m^, and 
M v i r based on angular-momentum conservation as outlined in MMW. In this standard model, the 
halo is assumed to contract adiabatically under the dissipative condensation of the gas following 
Blumenthal et al. (1986). 

A major novel ingredient considered here is the possible avoidance of the adiabatic contraction, 
hereafter AC (see § 4 for possible physical motivation). Another straightforward novelty adopted 
here is splitting the disc into a gaseous and a stellar component. Using Toomre's critical surface 
density Z cr i t , we compute the stellar mass as the disc mass with surface density £ > Z cr it, mimicking 
a star-formation threshold density in agreement with the empirical findings of Kennicutt (1989). 
Although this is clearly an oversimplification, it helps explaining (a) why the gas mass fraction is 
higher in lower mass galaxies (Kannappan 2004), (b) why the scale length of the cold gas is larger 
than that of the stars, and (c) the truncation radii of stellar discs (e.g. van den Bosch 2001). 

The parameters c, X and m& are assumed to be independent random variables, each with a 
log-normal distribution in which the mean could vary systematically with M v j r . In order to pur- 
sue the comparison to observations, we construct a large sample of random realizations of model 
galaxies. For a given M v j r , the halo concentration c is drawn from a log-normal distribution, with 
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Table 1: Slopes of observed versus model scaling relations. The upper 2 rows refer to the observations, 
while the lower 5 rows to the model. V bs is the observed (inclination corrected) rotation velocity. V3.2 is the 
model circular velocity at 3.2R S . L\ is the observed (dust corrected) /-band luminosity. M s and Mi are the 
model stellar and total disc mass, respectively. R s and Ri are the corresponding disc scale lengths. 

c = 12(V vir /100)" - 33 and Gi nc = 0.32, based on the CDM-simulated haloes (Eke, Navarro, & 
Steinmetz 2001; Bullock et al. 2001a). The baryonic spin parameter X is initially drawn from an 
independent log-normal distribution, with X = 0.035 and G\ n x = 0.5, again based on the simulated 
properties of CDM haloes (Bullock et al. 2001b), but we are forced to allow lower scatter values 
(see below). We then draw nii from an independent log-normal distribution, but consider nii (M v i r ) 
and Gi nmd as free parameters. Once we compute the disc mass and radius, we derive the stellar 
surface density profile and evaluate its mass M s , exponential radius R s , and "observed" circular 
velocity V3.2 at 3.2/? s . 

For a fair comparison with the data, we attempt to reproduce similar sampling in the model 
realizations. We first sample with equal probability in logM v i r , and then select a subsample for 
which the distribution of the derived V3.2 resembles that of the data. The model scaling relations 
and scatter are obtained by linear regressions in log space of V on M and R on M. 

3. Comparison of Models and Observations 

3.1 Slopes 

We start our exploration with the simplest MMW model possible. We keep c, X, and ma fixed 
at 12, 0.035 and 0.05 respectively, and compute the (scatter-free) scaling relations VM and RM. 
The apparent agreement between the predicted slopes and the /-band scaling relations, listed in 
lines 3 and 1 of Table 1 respectively, is misleading. This comparison naively assumes that all the 
disc mass is transformed into stars and with a constant Tj. A more realistic comparison of or model 
is with the relations listed in the second line of Table 1 , where M s replaces L/ in the data. Clearly, 
the basic MMW model yields steeper slopes than observed. The model VM slope gets smaller 
when the proper scaling of c with halo mass is incorporated, and it becomes almost as shallow as 
observed when the star-formation threshold density is imposed. A simultaneous match of the RM 
slope is obtained when nii is taken to scale with halo mass, nii oc M^ r 16 , as qualitatively expected 
from stellar feedback effects and as inferred from rotation curves (Persic, Salucci & Stel 1996). 

3.2 Zero-points 

Fig. 1 demonstrates that the model which matches the slopes in §3.1 suffers from severe zero- 
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Figure 1: Scaling relations for the model which matches the slopes (last row of Table 1, open circles), versus 
the mean and 2a scatter of the observed relations (dashed lines). The effects of 2a variations in m&, c and A. 
are shown, where 0i n „, d = 0.23, a\ nc = 0.32 and a^x = 0.50. The arrows indicate the direction of increasing 
rnj, c, and X. 



point offsets compared to the observed relations. This could have been fixed by a ~ 0.4 dex de- 
crease in Tj, but realistic IMFs limit the correction to 0.15 dex, which we adopt in our models 
below. Also seen in Fig. 1 is the impact of varying c, X. and over the 2-c range of their distribu- 
tions. Clearly, the VM zero-point is mostly sensitive to c, while the RM zero-point is most sensitive 
to X and to a lesser extent to the other variables (the secondary dependencies are partly due to the 
effect of AC). 

The RM zero-point could have been matched by decreasing X or increasing m<i, but these 
changes would have left the VM offset unchanged or worsen it. Furthermore, such changes would 
worsen the conflict with the observed, uncorrected residuals (see § 3.4 below), and thus fail to pro- 
vide a simultaneous fit. The VM zero-point could have been matched with the acceptable 0.15 dex 
decrease in Tj plus a factor of 2 decrease in c, 1 but this requires a drastic change of cosmological 
parameters, and also worsens the problem with the observed, uncorrected residuals. 

The above is a reflection of a well known failure of current semi-analytical models of galaxy 
formation, being unable to simultaneously fit the luminosity function (and hence the average virial 
mass-to-light ratio) and the TF zero-point. A VM zero-point match could be obtained if V bs — Kir. 
but realistic halo concentrations in the standard ACDM cosmology with o& = 0.9 imply Vobs/VW — 
1.4, and standard AC actually boosts it up to Vobs/Kir — 1-6. 

We are led to the conclusion that the only way to obtain a simultaneous match to the slopes 
and zero-points of both the VM and RM relations, within the realistic range of parameter values in 
the ACDM scenario, and with a hope to match the uncorrected residuals (§3.4), is by weakening 
or completely removing the effects of AC on V3.2 (see Fig. 3). 

3.3 Scatter 

The observed scatters about the scaling relations are G\ n (V\Lj) = 0.14 and Oi n (R\Lj) = 0.35. 
The measurement errors are estimated to be CTi n ~ 0.09 for each of the observables V, R, and Lj. 
In addition, there is a contribution to the scatter in our evaluated M s due to scatter in the relation 

'introducing a central density core does not affect V3 2 significantly (e.g. Dutton et al. 2003). 
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between color and Tj, which we estimate to be Gi n r, ^ 0.23 based on Bell et al. (2003). The 
remainder of the observed scatter is assumed to originate from the intrinsic scatter in c, X, and m&. 

Table 2 illustrates how each of these five sources of scatter contributes to the total scatter in the 
two relations. Note that although the different sources of scatter are assumed to be uncorrected, 
the residuals in the derived quantities (Lj, R s , V3.2) are partly correlated. This explains why the total 
scatter is less than the quadratic sum of the scatter in each component. 

As evident from Fig. 1 and Table 2, the scatter in X completely dominates the RM scatter. The 
observed RM scatter places a robust upper limit of 0\ n x ~ 0.25, independent of the other properties 
of the structural correlations. This is roughly half the value found in cosmological simulations, and 
may indicate that disc formation only occurs in a subset of haloes (with a restricted range in X). 

In the VL scatter, the contributions of observational errors and scatter in Tj and in c are all 
comparable, while the contribution of 0\ n x = 0.25 is smaller, and that of is negligible. While it 
would have been possible for scatter in Tj or c alone to account for all the intrinsic VL scatter, each 
would have resulted in a strong residual correlation, in conflict with the observations (§3.4). We 
therefore conclude that both c and T/ have to contribute comparably to the VL scatter. Note that 
with a significantly higher <5\ a \, the residuals about the VL relation would have been correlated with 
X. Given that the radius, and thus the surface brightness, is a strong function of X, this would have 
implied a surface-brightness dependence in the VL relation, contrary to observations (Courteau & 
Rix 1999; hereafter CR99). 
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Table 2: sources of scatter in the size-luminosity and velocity-luminosity relations. The intrinsic sources of 
scatter are: 0\ nc = 0.32; G\ B \ = 0.25or0.50; Oi nmd = 0.23; Oi n T/ = 0.23, and Oi ner r = 0.09. 

3.4 Residual Correlation 

As noticed by CR99, the degree of correlation between the VL and RL residuals at a given L 
provides a powerful constraint. Courteau et al. (2003) measured a weak anti-correlation between 
these residuals, fit in log space by a line of slope S = Alog V/Alog/? = —0.07 ± 0.03, confirming 
the original result of CR99. Note that S = for completely uncorrected residuals. 

The slope S is connected to the relative contributions of the disc and halo to the potential, 
and hence to the circular velocity V2.2, measured at 2.2/? s . The disc contributes V d 2 isc oc M s /R s , 
namely S = —0.5. A fixed halo, not reacting to the baryonic contraction, contributes V h 2 alo oc R s , and 
thus (S = +0.5) for an NFW density cusp. One can therefore expect that a weak anti-correlation, as 
observed, would be obtained with comparable contributions of disc and halo. AC increases the halo 
contribution to V3.2 in a way that anti-correlates with the disc radius, thus lowering S toward small 
or negative values. In order to maintain a small negative S with AC, the contribution of the halo, 
therefore, has to be larger. Indeed, assuming AC, CR99 interpreted their uncorrected residuals as 
indicating a "sub-maximal" disc with Vdi SC /V2.2 — 0.55. Matching the low negative S observed with 
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Figure 2: Correlations between residuals of VL and RL at a given L for our fiducial model with no AC in a 
halo of V v ; r = 100. Top-left: scatter in X (green circles) and c (blue squares) only. Bottom-left: added scatter 
in (red triangles). Top-right: added density threshold for star formation. Bottom-right: added scatter in 
Y/ (magenta stars). The 4 panels refer to residuals at a given M<j, M s and M s respectively. The linear fit 
(Alog V on A log R) is shown (solid line) and its slope is marked. The mean and 2a scatter about the observed 
relation are shown (dashed lines). 



our models that include AC, requires very low values of m&/\, pushing the parameters outside their 
reasonable ranges. Furthermore, this is at the expense of increasing the zero-point offsets (§3.2), 
so this model does not provide a simultaneous fit to the residuals and the zero-points. 

When the AC is eliminated, our model obeys the zero-point constraints and reproduces the 
weakly anti-correlated residuals. Fig. 2 illustrates how the various sources of scatter affect the 
distribution of residuals AlogV and Alog/? at a given M. When we use the total disc mass and 
radius, with scatter only in X and c, we still find a strong anti-correlation, S = —0.22. Adding scatter 
in nid reduces it to S = —0.20. When considering stellar mass and radius, we obtain a significantly 
weaker S = —0.13; this is because £ cr i t links M s to X (van den Bosch 2000, Firmani & Avila-Reese 
2000). When adding the scatter in Tj, there is another significant reduction to S = —0.087, due to 
the fact that the Tj scatter by itself induces a positive correlation of S = +1 (bottom-right panel). 
Finally, when adding the (uncorrected) measurement errors we obtain the residual distribution 
shown in Fig. 3, with S = —0.074 — in good agreement with the observations. 

Fig. 3 shows the joint distribution of structural parameters in our successful model compared 
to the observed results. The model assumes no adiabatic contraction. It adopts the c distribution of 
CDM haloes. It assumes X = 0.035, as for haloes, but G ln ^ = 0.25, about half the width of the halo 
distribution. It has m d = 0.07(M vir /10 12 M Q ) a16 , with G\ nmi = 0.23. The translation to M s is done 
using a T/ that corresponds to a Kroupa IMF, with Of l = 0.23. A measurement error of G\ n = 0.09 
is assigned to each of the three observed quantities. 



4. Conclusion 

We used the slope, zero-point, scatter and residuals of the observed velocity-mass and radius- 
mass distributions for disc galaxies to place constraints on the basic CDM-based model of galaxy 
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Figure 3: Distribution of structure parameters for our fiducial model (no AC, Kroupa IMF) versus the 
observations. The top-right panel refers to the residuals. Colored symbols refer to central surface density 
[M pc- 2 ]: E < 200 (black triangles); 200 < E < 630 (blue squares); 630 < E < 2000 (green circles); 
and 2000 < Eo (red stars). Dashed lines refer to the mean and 2a scatter in the data. 



formation. By applying these constraints simultaneously, we obtain non-trivial conclusions. Our 
simple model assumes independent log-normal distributions for 3 random variables at each halo 
mass: halo concentration c, baryonic spin X, and disc mass fraction m^. Additional sources of 
scatter come from the stellar mass-to-light ratio Tj and the measurement errors. 

Contrary to previous studies, we take account of the fact that discs are not purely stellar, and 
that the stellar mass-to-light ratio is not constant. For our model-data comparison, we have con- 
verted the observed /-band luminosities to stellar masses, using the relation between color and T/. 
This immediately demonstrates that even the slopes of the scaling relations are not automatically 
reproduced from the basic halo virial relation M v - ir oc V v 3 ir oc /j3 We find that for a simultaneous 
match of the observed slopes, given the c distribution as predicted for CDM haloes, stars have to 
form only above a threshold density E C n t , and the stellar mass fraction has to increase with halo 
mass. This is consistent with the expected effects of supernova feedback, known to be more ef- 
ficient in suppressing star formation in haloes of lower masses (e.g., Dekel & Silk 1986; van den 
Bosch 2002; Dekel & Woo 2003). 

The RM scatter is dominated by the scatter in spin parameter X. The relatively small scatter 
observed requires robustly that G\ n \ ~ 0.25, about half the value predicted for typical CDM haloes. 
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This could be due to a selection effect in the data (e.g. against LSB or SO galaxies), or may indicate 
that discs form in a special subset of haloes with a narrower X distribution. Interestingly, haloes 
with a quiescent history, which are the natural hosts of discs, indeed tend to have a narrower X 
distribution (D'Onghia & Burkert 2004). Alternatively it may indicate that the baryons somehow 
developed a narrower X distribution than that of the dark matter haloes. The VM relation is much 
tighter because it is largely independent of the scatter in X. The intrinsic VM scatter is determined 
by the scatter in both c andTj, as each alone violates the constraint from the uncorrelated residuals. 

Matching the zero-points of the scaling relations, combined with the weak anti-correlations 
between the residuals, turns out to be the greatest challenge. It could have been achieved by drasti- 
cally lowering c relative to the predictions for typical ACDM haloes, or by assuming an unrealistic 
top-heavy IMF. In both cases the parameters have to be pushed away from their reasonable limits, 
and the zero-points are matched without reproducing the uncorrelated residuals. One could alter- 
natively achieve weakly anti-correlated residuals by assuming unrealistically low m&jX, but only 
at the expense of worsening the zero-point offsets. 

The only way we found to match simultaneously the zero-points and the uncorrelated residuals 
is by eliminating the adiabatic response of the dark matter to the condensation of the gas into the 
central disc. Once the AC is removed, it becomes relatively straightforward to match simultane- 
ously all the constraints (though with a low scatter in X). Except for a model with unrealistically 
low stellar mass-to-light ratios, no other reasonable fix of the model parameters comes even close 
to simultaneously matching the zero-points and the residual constraints, making our conclusion 
quite robust. Our successful model with no AC predicts that discs contribute about 50% to V2.2 in 
the smaller, low surface brightness (LSB) galaxies, but as much as 80% in the brighter, high sur- 
face brightness (HSB) galaxies. Thus, whereas the LSBs are indeed sub-maximal, with the haloes 
dominating the inner regions, the HSBs are closer to maximum discs. 

In an ongoing work, we attempt to quantify the degree of contraction permitted by the data; 
yet to be confirmed is our impression that the corrected AC suggested by Gnedin et al. (2004) is 
not enough. Once the discs are not built by smooth spherical infall, the adiabatic contraction could, 
in principle, be counter-balanced by a variety of processes, and may even end up as an expansion. 
One such process is dynamical friction (e.g. El-Zant et al. 2004; Ma & Boylan-Kolchin 2004), 
which could be important if the discs are built by relatively big clumps (e.g. Mailer & Bullock 
2004). Indeed, the physics of gas infall in dark haloes suggests that most big discs have formed by 
clumpy, cold streams rather than smooth infall (Birnboim & Dekel 2003; Dekel & Birnboim 2005). 

Two caveats are worth mentioning. First, we should verify using more accurate color infor- 
mation that there is no systematic uncertainty in the conversion of Lj to M s , which could cause an 
overestimate of T/ for a given IMF. Second, our analysis assumed independence of the scatter in 
the variables at a given M. We suspect that c and X may be slightly correlated, e.g. via the halo 
accretion history (Wechsler et al. 2002), and one could imagine scenarios where ma is also slightly 
correlated with the others. To what extent these weak correlations may affect our conclusions 
remains to be seen. 

We thank A. Mailer & S.M. Faber for stimulating discussions. AD has been partly supported 
by ISF 213/02, NASA ATP NAG5-8218, a Miller Professorship at UC Berkeley, and a Blaise Pascal 
International Chair in Paris. SC recognizes the support of NSERC through a Discovery grant. 
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